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Abstract. Let (M, g) be an (n+l)-dimensional space time, with bounded curvature, with respect 
to a bounded framing. If (M, g) is vacuum, or satisfies a mild condition on the stress-energy tensor, 
", then we show that (M, g) locally admits coordinate systems in which the Lorentz metric g is well- 

' controlled in the (space-time) Sobolev space L 2 ' p , for any p < oo. 
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1. Introduction 



A well-known issue in the geometry of space-times is to understand the regularity of metrics 
with given bounds on the curvature tensor. This issue arises frequently in discussions and analysis 
of the behavior at the boundary and definitions of singularities for space-times, c.f. [11], [6], [7], 
[14] for example. 

More specifically, it has been an open problem for some time, cf. [4]- [6], [14] for instance, whether 
a space-time (M, g) which has curvature bounded in L°° in a suitable sense has coordinate charts 
| in which the metric g = g Q/ g is C 1 ' 7 n L 2,p , for any 7 < 1, p < 00. Here C fc ' 7 is the Holder space of 

functions whose k th derivatives are Holder continuous of order 7, while L k,p is the Sobolev space of 
functions with k weak derivatives in IP . 

The purpose of this paper is to provide an affirmative solution to this problem, at least for 
vacuum space-times or space-times satisfying a mild condition on the stress-energy tensor. 



The solution of the corresponding problem in Riemannian geometry has been known for some 
time, and it is useful to state the exact result in this context before considering the Lorentzian 
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(3JT), analogue. Thus, let (M,g) be a Riemannian n-manifold, with say C°° smooth metric g. Suppose 

there exists a point p E M such that 

(1.1) distg( P ,dM)>l. 

Let R = Rijki denote the Riemann curvature tensor of (M,g), and let B p (r) denote the geodesic 
ball of radius r about p in (M, g). Suppose one has bounds 

(1-2) \R\l°°(b p (i)) < C, vol g B p (\) > v , 

for arbitrary constants C < 00, v Q > 0. Then there exists a constant r a > 0, depending only on C, 
v a and n, such that the ball B p (r ) admits a coordinate chart U = {itk}, in which the metric g is 
C 1 ' 7 (~l L 2,p , for any 7 < 1, p < 00. Further, there exists a constant R , depending only on C, v , n 
and p, such that 

(1.3) \\9ij\\v* < Ro, 

where the norm is taken over the ball B p (r ). A proof of this result may be found in [12] for instance. 
By Sobolev embedding C 1 ' 7 C L 2,p , for 7 = 1 — jj*, so that (1.3) also gives a bound on gij in C 1,7 . 

A direct analogue of this result in Lorentzian geometry is false, due to the existence of large 
families of non-flat space-times for which the curvature norm |R| 2 = Rjj&zR^ vanishes identically. 
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Thus, consider for instance the class of vacuum plane- fronted gravitational waves on M 4 , with metric 
of the form 

(1.4) g = —dudv — h(x, y, u)du 2 + (dx 2 + dy 2 ), 

(1.5) *{x,y)h = 0. 

For such metrics, the two possible scalar invariants in the curvature tensor, namely 

|R| 2 = (R,R) = Rij-wRtf*', and (R, *R) = R ijkl (*R)^ kl , 

vanish identically. The vacuum Einstein equations impose only the condition (1.5), i.e. that h is 
harmonic as a function of (x, y). Thus, the function h may be an arbitrary function of u, and so is 
not controlled in any Holder or Sobolev space. It is thus clear that there is no coordinate system 
in which a general metric g of the form (1.4) is controlled in L ,p , or even 

To deal with this situation, one imposes bounds on the components of R in a fixed coordinate 
system or framing. An efficient way to do this is to choose a future-directed unit time-like vector 
T = eo and extend it to an orthonormal frame e a , < a < n, where the space-time dimension is 
n + 1. Since the space T 1 - orthogonal to T is space- like and 0(n) is compact, the particular choice 
of framing for T 1 - is unimportant. The norm of R w.r.t. T is then defined as 

(1-6) |R# = X^iw) 2 . 

where the components are w.r.t. the framing e a . 

Observe that if, at a point p € M, the vector T = T p is contained in a compact subset W of the 
future interior null cone Tj+M, then the norms (1.6) are all equivalent, with constant depending 
only on W. Hence, if K is a compact subset of the space-time (M, g) and T is a continuous vector 
field on K, then T, (or more precisely, Im T, where T is viewed as a section of the tangent bundle), 
lies within a compact subset of T + M, where T + M is the bundle of future interior null cones in the 
tangent bundle TM. 

To state the main result, we need the following definition, which is essentially just a normalization 
on the size of the region to be considered in (M, g), as is (1.1). 

Definition 1.1. Let Q, be a domain in a smooth Lorentz manifold (M, g). Then Q satisfies 
the size conditions if the following holds: The domain Q admits a smooth time function t, with 
Cq 1 < ||V£|| < c , for an arbitrary but fixed constant c Q < oo. Further, one has 

(1.7) d = B p (l) x [-1, 1] CC fi, 

i.e. the 1-cylinder C\ has compact closure strictly contained in f2. Here B p (r) is the geodesic r-ball 
about a point p in S, where S = So = i _1 (0) and the metric g on S is that induced from g. The 
product B p (l) x [—1, 1] is identified with a subset of O by the flow of Vi, i.e. (q, s) — > 7 ? (s), where 
j q (s) is the flow line of Vt, starting at q and terminating on the s-level set S s = i _1 (s). 
Let T = Vt/\\Vt\\ be the corresponding future-directed unit time-like vector field, and set 

(1.8) D = TmT\ Cl CC T+Cl 

The size conditions represent a Lorentzian analogue of the condition (1.1). They can always be 
realized by choosing f2 to be a sufficiently small open set in (M, g) and rescaling the metric up 
sufficiently. Essentially, they just serve to normalize the data. 

The main result of the paper is then the following: 

Theorem 1.2. Let Q be a domain in a vacuum (n + 1)- dimensional space-time (M, g), n > 2, 
satisfying the size conditions. Suppose that there are constants C < oo and v > such that 

(1.9) |R|r < C, vol g B p {\) > v Q . 
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Then there exists a constant r Q > 0, depending only on C, v Q , c Q , D, (and n), and a coordinate 
system (r, Xj), 1 < i < n, on the r Q -cylinder 

(1.10) C ro = D p (r ) x [-r ,r ] C C\, 

such that the components of the metric g Q/ g are in C 1 ' 7 fl L 2 ' p , for any 7 < 1, p < 00. Here D p (r) 
is the geodesic r-ball about p in the level set r = and the product structure is that induced by the 
flow of Vr. 

Further, there exists a constant R Q < 00, depending only on C, v Q , c a , D and the exponent p, 
such that, on C ro , 

(1.11) \\g a p\\ L 2, P < R . 
More precisely, for any k < 2, and < a, (3 < n, 

(1-12) \\d*g a/3 \\ Ll - k , P < R , 

where denotes any k-fold space-time partial derivative and the spatial Lx~ k,p norm is taken over 
any spatial slice {r = const} in C To . The constant R Q is independent of t in [— r ,r ]. 

The coordinates in Theorem 1.2 are geometrically natural; the time coordinate r is a Gaussian 
(equidistant) coordinate, while the spatial coordinates are chosen to be harmonic on the spatial 
slices {r = const}. It is not clear if there exist space-time harmonic, i.e wave, coordinates, in which 
g has this degree of regularity. 

The condition that (M, g) is vacuum, i.e. Ric g = 0, is used in a rather minor way. It is only 
used, via the Bianchi identity, to obtain LP bounds on the 2 nd time derivatives <9 T <9 T go Q of the 
components go a . This is equivalent to bounds on the 2 nd time derivatives on the components of 
the shift vector of the coordinates, i.e. a bound on the acceleration of the shift. Such acceleration 
components do not appear in any component of the curvature tensor R. All other estimates on 
dyga/3 are independent of the Einstein equations. 

The vacuum condition can be weakened to an assumption on the stress-energy tensor T in the 
Einstein equations, of the form 

(1.13) ||v T r|| L -i,p<c. 

Here L x 1,p is the dual space of Lo ,q , the space of Lx ,q functions of compact support on spatial slices 
{r = const} within C To , p^ 1 + q^ 1 = 1. 

The condition (1.13) will be satisfied automatically for many physically natural matter fields, cf. 
Remark 2.1. 

Theorem 1.2 is formulated in such a way that it is easy to pass to limits. Thus, suppose (Mj, gj) 
is a sequence of smooth space-times satisfying the hypotheses of the Theorem. There exist then 
domains J7j C (Mj,gj), points pi £ Jlj such that the size conditions (1.7)-(1.8) hold, with Di 
uniformly compact in T + Oj, (i.e. T does not become arbitrarily close to null cones). If (1.9) and 
(1.13) hold uniformly on fij, then there is a subsequence which converges to a limit C 1 ' 7 D L 2 > p 
space-time (M, g), defined at least on an r D -cylinder C To . Further, the convergence to the limit is 
C 1 ' 7 and weak L 2 ' p , and the bound (1.12) holds on the limit. 

Define a Lorentz manifold (M, g) to be weakly regular if g is a continuous Lorentz metric, with 

1 2 

g E L,' (M). It is well-known, cf. [11], [9] for example, that such metrics have a well-defined 
curvature tensor R in the sense of distributions. This leads to the following corollary. 

Corollary 1.3. Let (M., g) be a weakly regular Lorentz manifold, and let CC M be a domain 
with compact closure in M. Suppose the size conditions hold locally on Q, in that the constant \ 
is replaced by a small constant S so that B p (S ) CC M, for any pgfl. Suppose also the bounds 
(1.9) hold locally and uniformly on 5 Q cylinders as in (1.7), centered at any p £ fl. 
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Then 0, may be covered by a finite atlas of charts in which the metric g = g a p satisfies all the 
bounds in (1.12), except for the TP X bound on d T d T go a - The bounds in (1-12) depend, near d£l, on 
the distance of dQ, to <9M. 

If in addition the bound (1-13) holds distributionally on M, (e.g. (M, g) is a weak solution of 
the vacuum equations), then all bounds in (1.12) hold locally on (fi, g). 

We refer to the proof of Corollary 1.3 below for the precise meaning that (1.13) holds distribu- 
tionally. 

The proof of Theorem 1.2 and Corollary 1.3 follow in §2, while §3 concludes the paper with 
several remarks and extensions of these results, together with some open problems. 

2. Proofs of the Results. 

In this section, we prove Theorem 1.2 and Corollary 1.3. For clarity, the proof of Theorem 1.2 
is divided into several steps, each treating basically separate issues. In the following, as already 
above, space-time quantities are generally denoted in boldface while spatial quantities are not in 
boldface. 

Step I. (Initial Choice of Domain). 

Let B p (r) be the intrinsic geodesic ball about p in S. Since S is achronal in the cylinder C\ from 
(1.7), the extrinsic radius of B p (r) is bounded below for r small. Thus, if 7 is any space-like curve 
in (M, g) from p to x € dB p (ri) C S, for n small, then the length £(7) satisfies L{j) > l r±; the 
constants l Q and 77 depend only on c and D in Definition 1.1. 

Let V ri be the domain of dependence of B p {r\) in the manifold (M, g). Thus, by choosing 77 
sufficiently small, again depending only on c Q , D, one has 

(2.1) V ri CC n, 

i.e. T> ri has compact closure in $7. 

The region T> ri is globally hyperbolic and hence any pair of time-related points in T> ri may be 
joined by a time-like maximizing geodesic in T> ri . Recall from Definition 1.1 that the curves j x are 
the flow lines of Vi through x. For r 2 > small, (to be determined below), let 

(2-2) q = 7 P (-r 2 ), 

so q << p, i.e. q is to the past of p. For x to the future of p, x » p, let 

(2.3) t(x) = dist s (x, q) — dist s (p, q), 

so that r(p) = and t(x) > 0, for x » p. The distance r(x) is realized by a maximizing time-like 
geodesic o v (t) = exp q (r + t q )v from q to x; here v € T+M, with g(v,v) = —1, and r D = dist g (p,q). 
This normalization gives p = cr„ o (0), for some v Q <G T+M. Let 

N = Vr 

be the corresponding unit time-like vector field, so that N is the tangent vector to geodesies a 
issueing from q. Although N is well-defined and smooth along the individual geodesies <t„(t), (for 
any r until one reaches the boundary of (M, g)), at points where the exponential map exp q has 
cut or conjugate points, N is not uniquely defined. Of course, past such cut or conjugate points, 
the geodesies a v (r) are no longer maximal. Thus, in general, r is merely Lipschitz and N, as a 
vector field, is defined only almost everywhere to the future of q. At the end of Step I, (via the 
work in Step V), it will be seen that in fact r and N are smooth, in suitable domains of a definite 
size. In the following, unless stated otherwise, all geodesies are assumed to be maximal, i.e. they 
are not continued past the conjugate or cut points of q. 
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Next, let S = So = r l (0) and similarly let £ r be the r-level set of r in V ri . Since the geodesies 
are maximal, <j v (t) G S t , and p € S = So- Again, in general, S T is only Lipschitz. For r2 < n, 
consider the intrinsic geodesic ball # p (r 2 ) C S 1 and let 

(2.4) C = {x £ D ri : x = <r(r), r < r 2 , and <r(r) D £ p (r 2 ) / 0} C V ri . 

This is the "cone" of maximal geodesies a starting at q, hitting S within i? p (r 2 ), and terminating 
at time r = r 2 . 

Observe that the vector field ./V restricted to C stays within a compact subset of T + M. In fact, 
since N is parallel along its geodesic flow lines, this needs to be verified only at the base point q, 
where it holds by construction. It then follows from the curvature bound (1.9) and the remarks 
following (1.6) that 

(2.5) \R\ N <Ci = C 1 (C,D). 

Now the curvature bound (2.5) and the Rauch comparison theorem, cf. [2] for instance, imply 
that if r 2 is sufficiently small, depending (explicitly) only on C±, then the exponential map exp q 
restricted to the interior future null cone in T+, has no conjugate points in C. Thus, exp q is 
of maximal rank, and so a local diffeomorphism on C. In fact, for n is sufficiently small, again 
depending only on C and c G , no time-like or null geodesic within T> ri has conjugate points, and so 
exp x is of maximal rank on timedike geodesies in T> ri , Vx £ T> ri . 

Since V ri is globally hyperbolic and without timedike conjugate points, it follows then from [2, 
Thm.11.16] for instance, that any pair of points x,y £ V ri with y >> x, may be joined by a unique 
maximizing timedike geodesic in T> ri , provided T> ri is future 1-connected, (i.e. any pair of timedike 
curves joining x and y are homotopic through timedike curves). 

In general, T> ri need not be future 1-connected. Consider for example the past null cone of 
2-dimensional Minkowski space in hyperbolic coordinates t G ( — oo,0), <p € ( — 00,00), 

-dt 2 + t 2 d<j) 2 . 

If (f> is identified periodically, with any period, then the resulting space-time is globally hyperbolic 
but not future 1-connected. The future exponential map based at any point q has cut points; if the 
period of (f> is sufficiently small, or if q is sufficiently close to {0}, then cut points occur arbitrarily 
close to q. 

However, T> ri is future 1-connected if it has a simply connected Cauchy surface S, i.e. S = S ri — 
S (~l V ri , for S as in (1.7ff). To see this, let 71, 7 2 be two time-like curves with common endpoints 
in T> ri . The flow of the time-like vector field Vi gives a strong deformation retraction of 71 U 7 2 
onto a closed loop Ai U A 2 in S. If S is simply connected, then Ai may be deformed into A 2 within 
S. These two homotopies, time-like along Vi and space-like along the Cauchy surfaces St may 
be performed simultaneously, but with the latter at a larger speed than the former, to produce a 
time-like homotopy from 71 to 72. 

We will prove later that S, (or more precisely, a domain in 5 of a definite size), is simply 
connected. However, in order not to overburden the arguments to follow with such further issues, 
we assume in the following, through Step IV, that the Cauchy surface S C V ri is simply connected. 
This hypothesis will be removed in Step V, using the results obtained in the previous steps. 

It follows then that the exponential map exp q is a diffeomorphism onto C, when restricted to a 
suitable domain in T^M. The time function r is smooth in C \ {q}, as are the level sets S T n C, 
and there is a unique maximizing geodesic from q to any point in C. 

From now on, we consider S T C C, and so let S T denote the prior S T intersected with C. The 
level sets S T form a foliation of C by equidistant space-like hypersurfaces, with unit normal N. 

Step II. (Initial Curvature and Volume Estimates). 
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The geodesic congruence a on (C, g) satisfies the Riccati or transport equation 

(2.6) K' + K 2 + K N = 0. 

Here K = D 2 t is the 2 nd fundamental form or extrinsic curvature of the leaves E r , Rat is the 
symmetric bilinear form given by Rjv(-X") = (H.(N, X)X, N) and ' is the covariant derivative in the 
direction N. Hence, the bound (2.5) gives 

\K' + K 2 \ < C ± . 

This estimate holds on E T , for all r G [— r 2 ,r 2 ]. It then follows by standard comparison theory for 
the Riccati ODE (2.6) that if r 2 is sufficiently small, depending only on c D , D, and C\, then 

(2.7) \K\ L oo < C 2 , \K'\ L oa < C 2 , 

on all E r , r G [— ^,^2]. The constant C 2 depends only on r 2 and C\. The Gauss equation relating 
the curvature R of the ambient manifold (M, g) with that of the spatial slices E T reads 

(2-8) R-jjfcj = Rijki + K ik Kji - KuKj k , 

for spatial components (ijkl). This, together with the bounds (2.5) and (2.7), thus gives the bound 

(2.9) \R g \ L oo < C 3 , 

on E r , r G [— ^, r2]. Let cLK" be the exterior derivative, (w.r.t. the connection induced by <?), of K, 
when K is viewed as a 1-form with values in T£ r , i.e. dK(X, Y, Z) = (V X K)(Y, Z)-(V Y K)(X, Z). 
The Gauss-Codazzi equations are 

dK = R^, 

where K N (X,Y,Z) = (R(N,X),Y,Z). Hence, (2.5) also implies 

(2.10) \dK\ L °o < C 3 . 
We record also the well-known constraint equations: 

(2.11) SK = —dH — Ric(iV), 

R - \K\ 2 + H 2 = 2Ric(iV, N) + R, 
where H = tr K is the mean curvature and the operators <5 and d are taken on E T . 

Next, we use the bounds above to obtain a lower volume bound on the spatial slices S r , r € 
[— , 7*2], from that on the slice S = Sq in (1.9). To do this, let S = S n C. The domain 5 may be 
written as a graph over S = So C C via the time coordinate r in the usual way. Thus, each geodesic 
a = a v intersects S and 5 in exactly two points cr(ri), cr(r 2 ), with |rj| < r 2 . For x = o v {ji) G S, 
let n(x) = T2 — Ti, so that <r 1 ,(r 2 ) G S 1 . This gives a diffeomorphism <^> : S — > S, and hence 



wo/5= fdV§= I <f>*(dV d )= f JdVx, 



where J = detDcft is the Jacobian of 0. Since both S and E are space-like, the function u is a 
Lipschitz function, (cf. [11]), whose (weak) derivative is uniformly bounded, since both normal 
vectors ./V and T lie in compact subsets of interior null cones. In addition, the (uniform) time r 
exponential map, mapping E to E r has Jacobian uniformly bounded above and below on [— r 2 \, 
by the bound (2.7). (Recall that H = tiK measures the infinitesimal volume expansion or con- 
traction). It follows that the Jacobian J is uniformly bounded below, (depending only on C, D). 
Hence, 

(2.12) volT, > vi ■ volS. 
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Now the lower bound on volB p (^) C S in (1.9) does not immediately imply a lower bound on volS; 
(it could apriori happen that most all of the volume of B p (^) occurs outside S). However, in this 
case one can repeat all the estimates (2.6)-(2.9) when the construction of C is based at other center 
points q' in place of q. Thus, for p' G B p {^) C S, define q' as in (2.2) and let C be then as in 

(2.4). The same estimates as above then hold in C. The corresponding domains S' = SnC give a 
covering of B p {}-) C S. Hence the volume bound in (1.9) and the estimates above now do give the 
existence of points p Q £ B p (tj) C S such that 

(2.13) volZ Po >v 2 >0, 

where S Po is the level set of r, (i.e. r G ), containing p Q and v 2 = v 2 (v ,D,Ci). (The local estimate 
(2.12) does not in fact depend on the absence of future cut points of exp q , cf. the discussion 
concerning (2.59) below). 

Recall the standard volume comparison theorem in Riemannian geometry: if (N,g) is a Rie- 
mannian n-manifold, with Ric g > — (n — l)k, then the ratio 

volD(r) 
volD k (r) ' 

is monotone non-increasing in r. Here D(r) denotes the volume of a geodesic r-ball at any fixed 
point, while D^r) is the geodesic r-ball in the n-dimensional space form of constant curvature 
k. It then follows from the curvature bound (2.9), together with (2.13), that the geodesic balls 
D Po {r) C (E Po ,g) satisfy 

volD Po {r) > v 3 r n , 

for all r < r 2 , where v$ depends only on v 2 and C3. 

Observe that a similar estimate also holds for geodesic balls on other spatial slices S T , with 
So = S Po , for r € [— ^-,r 2 ]. Namely, the L°° bound on K in (2.7) bounds the infinitesimal 
distortion in the spatial metrics, and hence distances and volumes, under the flow of N. It follows 
that within the cylinder C Q centered at p a , the volume estimate above holds for balls D Pr (r) C S T ; 
thus, for p T = ct Vo (t), where p = ^(O), and for r < r 2 , one has 

(2.14) volD PT {r) > v 4 r n , 

= Vi(v , C, D). An upper bound on the volume of voID Pt (t) of the form (2.14) follows immedi- 
ately from the curvature bound (2.9). 

In the construction above, we have shifted the original base point p to a new base point p Q . 
However, one may now use these estimates to obtain equivalent volume bounds for the slices S r 
within the original cylinder C centered at p. This may be done by constructing a suitable chain, 
of bounded cardinality, of overlapping cylinders C\ from C Q to C. One then uses the arguments 
above on each Cj, together with the fact that upper and lower volume bounds of spatial slices are 
equivalent to upper and lower volume bounds of each cylinder Cj. 

Thus, in the following, we work on the original cylinder C from (2.4) centered at p; the bound 

(2.14) holds with p in place of p Q . 

Step III. (Local Coordinates). 

In this step, we define the cylinder C To and the local coordinate system on it, and obtain in 
addition some initial estimates on g a/ g. The local coordinates are Gaussian in time and harmonic 
in space, (Gaussian-harmonic coordinate system). 

Thus, the function r from (2.3) is chosen as the time coordinate on C. To construct spatial 
harmonic coordinates, start with the slice S = So within C. By (2.9) and (2.14), one has the 
bounds 

(2.15) \R g \ < C 3 , vol g D p (r 2 ) > v 4 . 
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It then follows, for instance from the discussion in §1, that there exists r Q > 0, depending only 
on C3 and V4, such that the geodesic ball D p (r ) C S admits a harmonic coordinate system {xi}, 
1 < i < n, in which the spatial metric g = g|s is controlled in L 2 ' p , i.e. 

(2.16) \\9ij\\i?* < Ro, 

where the L 2,p norm is taken on D p (r ), and R Q = R (Cs,V4,p). The harmonic functions X{ are 
solutions to the Dirichlet problem 

(2.17) AgXi = 0, Xi\ dD = 4>i, 

where D = D p (r ) and fa are suitably chosen boundary values, (approximating linear-type func- 
tions, cf. [12]). 

Let <pi tT = <j>i o ip Tj where Vv is the time r flow from S T to So along the integral curves of N. 
Thus, ip T maps a domain D T C S T diffeomorphically onto D and 4>i jT are functions defined on dD T . 
It follows that 

(2.18) N((f) i:T ) = 1 at dD T . 

Define the functions Xi on D T to be solutions to the Dirichlet problem 

(2.19) A gT Xi = 0, Xi\ dDT = (f> ijT . 

By (2.9) and (2.14), the estimate (2.15) holds uniformly on D T , for r <G [— ^,^2]. Hence, as 
with So, r D > may be chosen, depending only on C3 and V4, such that the functions {xi} form a 
harmonic coordinate system on D T C S T , on which one has the bounds 

(2.20) \\9ij\\v* <Ro, 

where the L 2,p norm is taken on D T and g = g T . The estimate (2.20) holds for all |r| < r a . 
This construction gives the local coordinate system (r, Xi), 1 < i < n, on the r -cylinder 

(2.21) C ro = D p (r )x[-r ,r ] 

about p, where the product structure is defined by the flow of Vr. For the remainder of the proof, 
S T is now redefined to be its intersection with C ro , i.e. S r = D T . 

The metric g in these coordinates has the form 

(2.22) g = (-1 + |£| 2 )(^) 2 + gijidxi + iidT)(d Xj + Zjdr), 
where ^ = {^} is the shift vector. Thus, 

d/dr = N + £, 

with N = Vr. The lapse function a of this foliation satisfies a = 1. 

On each slice S r C C ro , one has good spatial control, namely for gij = g^ = g|£ T , (2.20) holds. 
As usual, Latin indices denote spatial variables, i.e. 1 < i,j < n, while Greek indices a, {3 
denote space-time variables, < a, [3 < n. 

In the following, all Sobolev norms L k,p are understood to be spatial norms, i.e. the derivatives 
and norms are taken on spatial leaves S r . Thus, for emphasis or clarity, we sometimes write L^' p 
in place of L k,p . All estimates will be independent of r, for r < r . 

Step IV. (L 2 ' p Estimates of g a/3 ). 

In this next step, we extend the estimate (2.20) to include the remaining terms go a , < a < n, 
and also obtain estimates on the time derivatives of g a p. 

Before beginning, we first improve the estimate (2.7) on the 2 nd fundamental form. Recall the 
Simons', (or Bochner-Weitzenbock) formula, cf. [3, Ch. II]: 

D*DK = 5dK + d5K - K(K), 
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on (£ T ,<7 T ), where the term 7Z(K) is linear in the curvature and K; the exact form of 1Z(K) plays 
no role in the argument, but for completeness is given by 7Z(K) = Rico K + K oRic — 2RoK, where 
R o K is the action of the curvature tensor R on symmetric bilinear forms. The elliptic operator 
D*D = -trD 2 is the so-called rough Laplacian. 

In the following, we frequently write / G L k,p or / G Lx ,p as shorthand for / is uniformly bounded 
in L k,p along the spatial slices S T , |r| < r a . 

By (2.10), € and hence <5<ifr G Lx 1,p , for all p < oo; recall that these spaces are defined 
as following (1.13). Similarly, by (2.11), since dd = 0, 

(2.23) dSK = d(RiciV) G L~ 1>p ; 

here we recall that the operators 5 and d are spatial. The term 1Z{K) is also bounded in L°°. 
Hence, one has 

(2.24) D*DK = Qi, 

where Q\ is uniformly bounded in Lx~ 1,p , for any p < oo, while K is uniformly bounded in L°°. By 
(2.20), the coefficients of D*D in the local coordinates are controlled in L 2,p C C 1,7 . It then 
follows from standard elliptic regularity theory, cf. [10], that 

(2.25) \K\ T i, P <C 2 , 
where C2 = Ci(C, v a , c Q , D,p), on all spatial slices S T , \t\ < r Q . 

Spatial Estimates. 

Here, we prove that the components go Q also satisfy the L 2,p = L x ' p estimate (2.20) uniformly 
on E T . One has N = g 0a d a , Vx { = g ia d a . Hence, N( Xi ) = (N, Vxi) = g 0a g il3 g a p = g 0i . 

To obtain estimates on iV(xj), differentiate the harmonic coordinate condition (2.17), in the 
normal, (i.e. iV), direction. Let x[ = N(xi). Since Ax{ = 0, a standard computation, cf. [3, 1.184] 
for example, gives 

(2.26) Ax- = -A'xi = (D 2 Xi , K) - (dx t , SK + \dH). 

Here, as above and in the following, all metric quantities in (2.26) are on spatial slices S T . 

By (2.25), K is uniformly bounded in L 1,p . The term D 2 Xi is also uniformly bounded in L 1,p , 
since by (2.20) the spatial metric is uniformly bounded in L 2 ' p and hence the coordinate functions 
are uniformly bounded in L 3,p . Further, both SK and dH are uniformly bounded in L p . Thus, 

Ax'i = Q 2 , 

where Q2 is uniformly bounded in L p , |r| < r . As before, the coefficients of A are controlled in 
L 2 ' p . Further, by construction, cf. (2.18), x\ = 1 on <9£ r . Hence, standard elliptic regularity again 
gives 

(2-27) K\\v* = ||g°iU2, P <C 4 , 

where C4 = C^[C, v a , c Q , D,p). Observe also that 

(2.28) g 00 = -1, 

(since the lapse function a = 1). Hence g 0a G L 2,p , i.e. the L 2,p norm of g 0a is uniformly bounded, 
< a < n. 

From this and (2.20), it is then an elementary exercise in linear algebra to see that 

(2.29) \\ga(s\\ L ^ < C 5 . 

Briefly, g° 7 = (detgap)^ 1 Aq^ , where Aq^ is the (0,7) cofactor in the matrix g a p- The cofactor 
Aqq involves only gij, and hence by (2.20), Aqq G L 2,p . Thus detg a p G L 2 ' p . The same reasoning 
on g 0a then gives Ao a G L 2 ' p , for all a. Each determinant Aq^ may be expanded along the first 
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column to obtain a linear form in the variables goj, with coefficients (n — 1) x (n — 1) determinants. 
Thus, one has a linear system of n equations in n unknowns go«. The matrix of this system is the 
(n — l)-compound G n -\ of the matrix [gij], i.e. (G n -i)ki = detA^i, where A^i is the (k,l) cofactor 
of [gij]. Since [gij] is non-singular, and since [g^] non-singular implies that G n -\ is non-singular, 
(cf. [8, §1.4] for instance), it follows that this linear system is invertible. 

The components go« are rational expressions in {g^} and {A^}, each of which is now bounded 
in L 2 ' p . Hence goi is bounded in L 2,p . Finally, since 1 = g 0a g a o = g°°goo + g°*gi0> it follows from 
(2.28) that goo is also bounded in L 2,p . This establishes the bound (2.29). 

Recall that d T = N + £, while (N,di) = 0, i > 0, by construction. Since (di,d T ) = g«o is bounded 
in L 2 ' p , it follows that (£, di) is bounded in L 2 < p . Hence the shift vector £ = is bounded in L 2 ' p , 



(2.30) m\ L ^<c 6 . 

This completes the L 2 ' p estimates of g a/ g in spatial directions. 
1 st Time Derivatives. 

Next we turn to estimates on the time derivatives of g Q| g, i.e. L 1,p estimates for d T g a /3. To begin, 
using the Leibniz rule, and the fact that [d a , d T ] = 0, it suffices to estimate 

(2-31) (VaA.fy) = (Va a iV,^) + (V d J,d p ). 

Suppose first a > 0, (3 > 0, so (a, (3) = (i, j). The first term in (2.31) is then Kij, which is bounded 
in L 1,p by (2.25), while the second term is also bounded in L 1,p by (2.30), (and the L 2,p spatial 
bounds on g a p in (2.29)). This gives uniform L 1,p bounds on d T gij, i.e. 

(2.32) \\d T g l3 \\ L i, P <C 7 . 
It follows of course that also 

\\dx k d T gij\\ L p < C 7 . 

The bounds on <9 r go Q require more work. Writing d T = N + £ as above, the L 2,p spatial estimates 
above imply that £(go«) is bounded in L l,p , so one needs to obtain L l,p bounds on iV(goa). 

Recall that x- = iV(xj) = g 0i . Hence iV(g 0i ) = NN(xi) = x". To obtain estimates on x", 
differentiate the equation (2.17) in the N direction twice. This gives 

(2.33) AxJ' = -(2AV i + A"ar i ). 

Here, as before, all metric quantities are on the spatial slices S T . It has already been proved that 
x\ £ L 2 ' p . From the form of A' in (2.26), one then easily sees that 

A'4 G L p . 

Next, one has 

(2.34) A"xi = N(D 2 Xi, K) - N(dx t , 5K + \dH). 

To estimate these terms, let e a be a local orthonormal basis on S r , with V eQ e^ = at any fixed 
point in S r . Then the first term in (2.34) may be written 

N(D 2 Xi , K) = N(D 2 Xi(e a , e b )) ■ K(e a , e b ) + D 2 Xi (e a , e b ) ■ N(K(e a , e h )). 

By (2.25) and (2.7), K ab G L liP and NK ab G L p , while by (2.20), D 2 Xi G L 1,p . Thus 

(2.35) D 2 x t (e a , e b ) • N(K(e a , e b )) G L p . 
Further N(D 2 Xi(e a , e b )) = N(V ea dxi,e b ) = (V N V ea dxi,e b ), so that 

(2.36) N(D 2 Xi (e a , e b )) = {V ea V N dx h e b ) + {R(N, e a )d Xi , e b ) - (V eb dx h V £a N). 
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The curvature term in (2.36) is bounded in L°°, while the last term equals (D 2 Xi(et,), K(e a )), which 
is bounded in L°°. For the first term, one has 

(2.37) {V ea V N dxi, e b ) = e a {V N dxi, e b ) = e a {V eb dXi, N) = -e a (K(dxi, e b )). 

Since K € L 1,p , this term is bounded in LP . Combining these estimates, it follows that the first 
term in (2.34) is bounded in L p . 

For the next term in (2.34), N(dxi,dH) = (\7 N dxi,dH) + (dxi,V N dH) = {V dH dxi,N) + 
(N, V dXi dH), so that 

(2.38) N(dxt, dH) = -2K(dH, d Xi ), 
which is bounded in L p . 

Finally, -N{dxi,5K) = Ne a K{dx i ,e a )-NK{V ea e a ,dx i )-NK{e a ,V ea dx i ). By (2.7), the latter 
two terms are in L°°. For the first term, write Ne a K(dxi, e a ) = e a NK(dxi, e a ) — (V ea N)(K(dxi,e b )). 
The latter term here is —[K(e a )](K(dxi, e b )), which is bounded in L p by (2.25). For the first term, 
since NK(dxi,e a ) is bounded in L°°, e a NK(dxi,e a ) = div(N(K(dxi))) is bounded in L~ 1,p , since 
the derivatives e a are spatial. This shows that 

(2.39) N(dxi,5K) £ L~ l ' p . 

Thus, combining these estimates on (2.33) gives a uniform bound on Ax" in L~ l,p . On the 
boundary <9X T , one has x'( = 0. It then follows from elliptic regularity as before, (as in (2.24)), that 

N(g 0i ) = x'l G L 1 *. 

Of course, by (2.28), iV(g 00 ) = 0. As above, go« is a rational expression in {gij} and {^4ofc}- The 
bound (2.32) implies that each of these has iV-derivative in L l,p and hence, by the same arguments 
as before, -/V(go a ) G L 1,p . This gives uniform bounds 

(2.40) ||<9 rga/3 || L i,p < C 8 . 

This completes the estimates for the first time derivatives on spatial slices. In particular, all 
Christoffel symbols are bounded in L^ p . 

2 nd Time Derivatives. 

Finally, we obtain LP estimates on the 2 nd time derivatives d T d T g a/ 3. To do this, take d T of the 
term (Vd a d T ,dp) in (2.31). One then obtains 

(v 9a VaA, dp) + (R(d T , d a )d T , dp). 

The curvature term is bounded in L p , in fact L°°. (This uses the fact that £ is controlled, so the 
framing is controlled). 

Write Vg T d T = Fg 3 7 . Hence 

(v da v dT d T ,dp) = rl (v da d^dp) + d a (rl )(d J ,dp). 

By the first derivative estimates above, the Christoffel symbols are bounded in Ll' p . The first 
term is a product of Christoffel symbols, and hence is bounded in Lx' p ^ 2 C L n p/( 2n -p) ; by Sobolev 
embedding. For p > n, (recall that p is arbitrarily large), L n p/( 2n -p) ^ £ p . Thus the first term is 
bounded in L p . 

For the second term, if a > 0, then d a (TQ ) is bounded in L p by (2.40). Hence, this gives 

(2.41) \\d T d Tgij \\ LP < C 9 , 

uniformly in r. It remains to estimate the second time derivatives of go Q . These correspond to the 
second order time behavior of the shift vector £. 
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These estimates are the most involved, and are the only estimates dependent on the Einstein 
equations. To obtain these estimates, one needs to differentiate (2.17) three times in the normal N 
direction. Thus, from (2.17) again, 

(2.42) AsJ" = -3A'x'{ - 3A'X - A'"^. 

Recall x\ = g°* is bounded in L 2,p , while x'[ = A^(g *) is bounded in L 1,p . From previous work, it 
is then straightforward to bound the first two terms on the left in (2.42). 
To see this, one has 

A'x'l = {D 2 x'l, K) - {dx'l, 5K + \dH), 

which is bounded in L^ 1,p , since x'[ G Ll' p . To estimate the term A"x-, one just replaces Xi in the 
estimates (2.34)-(2.39) by N(x{). Using the fact that N(xi) <G L 2,p , one sees by checking term by 
term that this is bounded in L~ l > p . 
It remains to analyse 

(2.43) A"'xi = NN(D 2 Xi , K) - NN{d Xi , 5K + \dH). 
To begin, 

(2.44) NN(D 2 Xi, K) = (V N V N D 2 x u K) + 2{V N D 2 x h V N K) + (D 2 x h V N V N K). 

By the Riccati equation (2.6), Vat if is bounded in L°°, while by the estimate on (2.36), VjyD 2 Xi 
is bounded in L p . Hence, the middle term in (2.44) is bounded in L p . For the last term, taking the 
A^-derivative of the Riccati equation (2.6) gives 

(2.45) {D 2 x u V N V N K) = -2(D 2 x h (V N K)K) - (D 2 x u V N R N ). 

The first term in (2.45) is bounded in L°°. The second (curvature) term in (2.45) will be analysed 
below. 

For the first term in (2.44), one needs to take iV-derivatives of all terms following (2.35) to (2.37). 
The only one which is not bounded in L~ 1,p by previous estimates is the term 

(2.46) N(R(N,e a )dxi,e b ) = (V N K)(N,e a ,dxi,e b ) + lower order. 

We return again to this curvature term below, and proceed with the second term in (2.43). 

Again, to estimate this, one takes A^-derivatives of the estimates in (2.38) to (2.39). This gives 
first 

NK(dH,dxi) e L p , 

since NK is bounded. For the 5K term, modulo lower order terms, this is of the form 

N(V N V ea K(dxi), e a ) = N(V e y N K(d Xi ), e a ) + N(K(N, e a )K(d Xi ), e a ) = 

= (V N V ea V N K{dxi), e a ) + (VatR)(AT, e a , K(dxi),e a ) = 

= (V ea V N V N K{dxi), ea) + (VjvR)(AT, e a , K{d Xi ),e a ) + (R{N, e a )V N K(d Xi ), e a ); 

the equalities here are understood to be modulo lower order terms. Modulo terms bounded in L p , 
this may be rewritten as 

(2.47) V ea (V^R)(AT, K(dxi),N, e a ) + {V N R)(N, e a , K(d Xi ), e a ). 

Combining these estimates gives then uniform L x l,p bounds on all terms in (2.43), except for 
the four curvature terms of the form (VatR)A^. 

We obtain bounds on the curvature terms via the contracted 2 nd Bianchi identity on (M, g): 

dli = -dRic, 

or more precisely, cf. [3, 16.3], 

(2.48) 6R{X, Y, Z) = -dRic(y, Z, X). 
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Write, on M, 

(2.49) (V iV R)A r = -5R + (5R = dRic + (5R, 

where the divergences 3, 5 are the space-time and space-like divergences on E r respectively. 

Now the space-time curvature R is bounded in L°°, (in bounded framings). Hence the spatial 
divergence <5R is bounded in Lx ,p , for any p < oo. Thus, the curvature term (2.46) may be 
rewritten, modulo Lx 1,p , as 

(VjvR)(A?", e a , dxi, e b ) = dRic(dx«, e b , e a ) G L~ 1,p , 

where the last estimate follows since all derivatives of Ric are taken in spatial directions. Similarly, 
for the second curvature term in (2.47), one has, for the same reasons, modulo L x ' , 

{V N K)(N, e a , K(dxi), ea) = dRic(K{d Xi ), e b , e a ) G L' 1 ^. 

This leaves left the two curvature terms: 

(2.50) (D 2 Xi ,V N R N ) and V ect (V^R)(A^, K(d Xi ), N, e a ). 

The second term is of form 5d~Ric(N, K(dxi)), (to leading order), which cannot be controlled 
without the Einstein equation, since it involves differentiation in the N direction. Similarly, 

(2.51) VtvRjv = dRic(iV, •,•), modulo L~ 1,p , 

is not controlled without the Einstein equations. 

Dropping the usual constants, the Einstein equations on (M, g) are 

(2.52) Ric -*9 = T, 

where T is the stress-energy tensor. Suppose (M, g) is vacuum, or more generally, suppose the 
stress-energy tensor T satisfies 

(2.53) dT(N,e a ,e b ) £ L~ 1,p , 

where e a ,et, are spatial. Since the Einstein equations and the bound (2.5) imply that V ea T(N) € 
Lx~ 1,p , (2.53) is equivalent to 

(2.54) VjyT G L x ~ 1,p , or C N T G L~ 1,p , 

where Cn denotes the Lie derivative in the direction N. 

Combining the bound (2.53) with the estimates obtained above on the terms in (2.42) then gives 

(2.55) Ax'f = ANN(g 0i ) G L~ 2 ' p , 

since 5dRic(N,K(dXi)) G L~ 2 ' p , not Lx P . As before, since the coefficients of the Laplacian 
are well-controlled, (i.e. bounded in L 2,p ), and the functions x'[' have boundary values, elliptic 
regularity gives 

iViV(g 0i ) G L p . 

(One sees this from duality in the standard way, using the fact that A : L 2,p — > L p is an isomorphism, 
cf. also [13]). 

Applying the linear algebra argument as before then gives the bound 

(2.56) d T 8 T g 0a G L p . 

This completes all of the estimates on g a p. We refer to Remark 2.1 below for further discussion 
on use of the Einstein equations. Finally, the assumptions (2.53) or (2.54) are equivalent to (1.13). 
To see this, both T and N lie within a compact subset of T + f2. Since the covariant derivative w.r.t. 
T or N involves only the pointwise behavior of these vector fields, one may replace T by N in 
(1.13), which then corresponds to (2.54). 
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Step V. (Issue of Cut Points). 

In this final step, we show that by passing to a smaller cylinder if necessary, of definite size within 
C ro , the exponential map exp q > at a suitable base point q', has no future cut points. By the work 
above, this will complete the proof. 

To begin, return to the "cone" C in (2.4). Let V be the collection of time-like unit vectors in 
T q T> ri for which there is a maximal geodesic a v issueing from q and terminating in C. Let 

(2.57) C = {av:v€V,s<n}c T+V ri , 

so that in particular exp q C contains C. Here, as in Step I, r\ is chosen so that no geodesic cr v (s), 
for v € V, has conjugate points with s < r±. Of course the geodesies cr v (s) now are no longer 
necessarily maximal. Since exp q is of maximal rank on C, we work on the pullback (C,g), where 

(2.58) g = (ex Pq )*g. 

The domain C is a compact connected cone in T g M, w.r.t. the vector space structure. By the 
Gauss Lemma, the straightline generators of this cone are geodesies in the g metric. Let r denote 
the distance to the origin {0} w.r.t. g, within C. This is now a smooth function on C \ {0}, and 
serves as the parameter for the geodesies from {0}. Inside the cutlocus of exp q on C, (i.e. where 
exp q is a diffeomorphism) , r is just the lift of the function r from (2.3), up to an additive constant. 

The level sets c C of r are smooth, and hence the images exp q T,^ are smoothly immersed 
submanifolds in C. The original Lipschitz level surface S T is just the part of exp q T,^ contained in 
the domain U T = {x : r(x) > r} n C. Observe that the timelike exponential map exp of (C,g) 
based at is a diffeomorphism onto C; this map has no conjugate or cut points within C. 

Thus, we are now in exactly the same situation as at the end of Step I, with (C, g) in place of 
(C,g). Set S = (expq)' 1 ^) C C, so that S is an embedded hypersurface in C. Since exp q is a local 
isometry of S onto S, one has 

(2.59) VOlgS > VOlgS. 

Hence, by the same reasoning as in Step II, the volume estimate (2.14) holds on the smooth 
hypersurfaces Sy. 

We may thus apply the work in Steps II - IV to conclude that there is an r Q > 0, depending only 
on C, D, c , v , and an r -cylinder C To C C, centered on p, on which there are coordinates (r, Xj) 
in which the metric g is controlled in L 2,p in the sense that (1.12) holds. 

It remains to prove the existence of a suitable cylinder C' "downstairs", i.e. within (fi,g), with 
these properties. Thus, let now denote the part of T,^ contained in C To . Each is an n-ball, 
topologically. Recall, as in Step II, that S n C To is a graph over S^ o , for some r Q . Hence, S n C To 
is also an n-ball topologically. By construction, both and S n C To have a uniform lower bound 
on their volume and size, i.e. distg(p, d(S n C To )) > > 0, (and similarly for T,^ o ). These bounds 
depend only on the initial bounds on C, v Q and D, c . 

Let U be the interior of the cutlocus of exp q \^. The domain U is starshaped w.r.t. the origin 

{0} in T+P ri . Further, again by construction, S n C To is contained in the closure U of U. If 
(S n C To ) n dU / 0, one may perturb it slightly, along the geodesic straight lines to {0}, to obtain 
an n-ball S' C C To with S' C U. As before, the ball S' has a definite lower bound on its volume 
and its size. 

The exponential map exp q now gives a diffeomorphism, in fact an isometry, from (S", g) to (S 1 , g), 
S' = exp q S'. Let V be the domain of dependence of S' in (M, g). This gives a globally hyperbolic 
region V' C (M, g), with a simply connected Cauchy surface of definite size and volume. The work 
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of Steps I - IV may now be applied to this situation within (M, g) to produce a new cylinder C', 
centered at p, satisfying the bounds (1.12). 

This completes the proof of Theorem 1.2. ■ 

It is an open question whether Theorem 1.2 holds without an assumption of the form (1.13), i.e. 
whether there exist coordinate systems in which (1.12) holds under only the bounds (1.9). 

Remark 2.1. From the physical point of view, most stress-energy tensors T derive from matter 
fields satisfying a hyperbolic system of PDE, of 1 st or 2 nd order. In such a situation, these equations 
can frequently be used to interchange a time derivative VatT on spatial slices, with a spatial 
derivative VxT, modulo lower order terms, e.g. (VjvT)(X) ~ (^7xT)(N), modulo lower order 
terms. This is exactly the process used via the 2 nd Bianchi identity above. For instance, this is 
easily seen to be the case for electromagnetic fields, via use of the Maxwell equation dF = 0. 

When the matter equations allow for such time-space replacement, modulo lower order terms, 
the condition (1.13) is of course not necessary in Theorem 1.2. 

Next we turn to the proof of Corollary 1.3. 

Let (M, g) be a weakly regular space-time, satisfying the size conditions, and satisfying the 
bound (1.9) distributionally. Thus, the components of R, well-defined as distributions, are in fact 
bounded in L°° . 

Any such space-time (M, g) is a limit of a sequence of (C°°) smooth space-times (M, g k ), cf. [9, 
Theorem 4] for instance. The metrics g& are obtained in the usual way by taking the convolution 
of g with a sequence of smooth mollifiers. The local size conditions and local volume bound in 
(1.9), with d in place of ^, depend only on the C° behavior of the metric. Since the convergence to 
the limit is C°, it follows that local size conditions and local volume bounds hold uniformly on the 
sequence (M, g fc ). Similarly, the fact that |R|t is bounded on (M,g fc ) implies that the curvature 
|Rg Jr fc °f (M, g fe ) is uniformly bounded, for unit time-like vector fields — > T, as k — > oo. 

It then follows from Theorem 1.2 that for any p £ H, there are r G -cylinders ((C ro CC 
(M,g fc ), with ph — > p, and coordinates on {C To )k in which the metric g^ is controlled: this in the 
sense that the bounds (1.12) hold, with the exception of the bound on <9o<9ogOa- Since the bounds 
on C, v and D, c Q , hold uniformly on (M, g fc ), and p\. remains a bounded distance away from <9M, 
it follows that there is a limit cylinder (C ro ,p) C (M,g) on which (1.12) holds, again except for 
the bound on <9o<9ogoa- 

This proves the first part of Corollary 1.3. For the second part, one needs to make sense of the 
condition (1.13) on (M, g). It suffices to do this locally, i.e. on cylinders C To C (M, g). Of course 
VjvRic, or equivalently V^T, is well-defined in L~ l,p (Q), (since Ric <G L°°), but one needs to 
define it in L- 1 ' p (S r ). 

To do this, from the above, we know that g T = g|s T is in Lz P (£ T ), for any S T C C To . Let h be a 
symmetric bilinear form in L 1,q (Y, T ), of compact support, for a given r. Extend h into C To by the 
flow of N, so that C^h = 0. Hence, h is defined on all S r C C To . Formally, or alternately on the 
smooth approximations (M, g fe ) of (M, g), one has 

(2.60) f (h,V N mc)dV = I N{h, Hic)dV — f (V N h, Hic)dV. 

Observe that V N h = C N h + h(V 9a N,d b ) + h(d a ,V db N) = h(K(d a ),d b ) + h(K(d b ),d a ). Since 
K £ L 1,p (E r ) and p is large, it follows V xh is well-defined in L q (T, T ), for any r. Hence, the second 
term in (2.60) is well-defined, (since Ric is bounded in L°°). 
For the first term on the right in (2.60), one has 

(2.61) [ N{h, Kic)dV = / (h, Ric)dV - [ {h,Bic)HdV. 
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The second term in (2.61) is well-defined, since again H G L 1,p C C a . Thus, to define (1.13) on 
C To , we require that the derivative 

d f d f 

(2.62) — / (h,Kic)dV, or equivalently, — / (h,T)dV 

dr Jy, t dr J^ t 

exist for all r G [— r ,r ], for any h as above. Under this condition, the bound (1.13) on C To is then 
equivalent to a uniform bound on (2.62), for all r G [— r a ,r ]. 

Given this definition of the condition (1.13) on the limit C ro , it follows from the same proof as 
in Theorem 1.2 that <9o<9ogOa * s uniformly bounded in L P (S T ), |r| < r Q . This completes the proof 
of Corollary 1.3. ■ 



3. Concluding Remarks. 

We conclude the paper with several remarks extending the validity of Theorem 1.2, together 
with some open problems. 

Theorem 1.2 gives the existence of cylinders C of a definite size in the interior of M, and 
a definite distance away from any boundary 3M, on which there exist coordinates in which the 
metric is controlled in L 2,p . Of course, by rescaling up suitably to realize the size conditions, applying 
Theorem 1.2, and then rescaling back down, the cylinders C may be chosen to be arbitrarily close 
to dM.. (This is already implicit in Corollary 1.3). However, the coordinates may change in the 
(smaller and smaller) cylinders as one approaches <9M. 

The main reason it is necessary to stay a definite distance away from <9M in the proof is that 
one needs to work in globally hyperbolic regions, as in (2.1), which are future 1-connected. If M is 
globally hyperbolic and future 1-connected to begin with, then it is no longer necessary to stay a 
given distance away from 9M. We describe an alternate version of Theorem 1.2 in this context. 

Thus, suppose (M, g) is globally hyperbolic and future 1-connected; (for instance, (M, g) has a 
simply connected Cauchy surface). Choose any point q G M and time-like unit vector N Q G T+M. 
Let V C T+M be the compact cone of unit vectors forming a fixed angle 6 with N a , and let 
C = expgTV, v G V, be the corresponding "cone" in (M, g). Also, let N be the corresponding 
vector field tangent to the geodesies a v (r) = exp q (rv), and let S T be the level sets of r in C. We 
make the following assumptions: there exist constants C < oo and v > such that 

(3.1) |R|jv < C, within C, 



(3.2) volT, T > v ■ (diamT lT ) n . 

Although (3.2) is understood to hold for all r < r\, for a fixed r±, the methods used in Step II 
show that it suffices to assume (3.2) holds for a fixed r G > 0, where the size r depends only on C 
in (3.1). The bound (3.2) is of course scale-invariant. 

We then have the following result, valid at least up to a point in <9M. 

Corollary 3.1. Let (WL, g) be a globally hyperbolic, future 1-connected space-time satisfying the 
assumptions (3.1)-(3.2). Suppose also (M., g) is vacuum, or more generally, satisfies (1.13). Then 
there exist small constants d Q > 0, and r Q > 0, depending only on C, v Q , N Q and 6, such that if 

(j No {t) G M, Vr < d , 

then there exist coordinates (r, x«) on the cylinder 

(3.3) C ro {di) = D p {r d ) x [f C M, 

in which the metric satisfies the bounds (1.12). Here d\ is the largest value such that d\ < d Q and 
C ro (di) is contained in M. 
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The proof of this result is exactly the same as that of Theorem 1.2. In fact it is simpler, since 
the issues of global hyperbolicity and future 1-connected are assumed, and one works with iV in 
place of the vector field T from Definition 1.1. Thus, for d Q small as above, and for C(V) the cone 
on V in T g M, the exponential map exp q restricted to (C(V) \ {0}) n Bo(d\) is a diffeomorphism 
onto (C fi B q (dij) \ q, where B q {d\) = expq(Bo(di)) . The domain C fl B q (d\) plays exactly the same 
role as C in Step I. All the estimates of Steps II-IV then proceed just as before and are uniform for 
r€[^M u 

Remark 3.2. One may also derive a version of Theorem 1.2, (or Corollary 3.1), without the lower 
volume bound in (1.9). Thus, one has a uniform spatial curvature bound (2.9), either on the slices 
S r C C downstairs, or on the slices C C upstairs in T g M. (The bound (2.9) does not require any 
apriori volume bound). If the volume S T or is very small, then the injectivity radius of g T is very 
small, i.e. the spatial metric g T is highly collapsed in the sense of Cheeger-Gromov. However, (as 
before), the Rauch comparison theorem implies that the intrinsic exponential map on S T or is 
still of maximal rank on geodesic balls of a definite size, (depending only on the curvature bound). 
Suppose first one works in the situation of S r C C. Then just as before in time-like directions, 
one can lift the metric g T up to the tangent space by pulling back by its exponential map, i.e. 
consider the metric g T = (exp PT )* g T defined on balls in T Pt T, t . The metric g T has a uniform lower 
bound on the volumes of small balls. All the arguments in Steps II-IV can then be carried out as 
before on this "unwrapped" space-time and the corresponding unwrapped cylinder C To . This gives 
coordinates (r, Xi) for C To on which the metric g satisfies the bounds (1.12). The same procedure 
holds when working with Ey. These coordinates upstairs within T 9 M give then "multi- valued" 
coordinates downstairs in (M, g). 

To conclude, we mention two open problems. First, for certain purposes, the L°° bound on the 
curvature in (1.11) may be viewed as too strong. It would be of interest to know if a version of 
Theorem 1.2 holds with suitable L p bounds on |R|t> V > n /2, in place of L°° bounds. 

Second, it would be very interesting if the Einstein equations could be used to remove the 
dependence of these results on bounds on the full curvature, i.e. if bounds on the full curvature 
could be replaced by bounds on the Ricci curvature, possibly introducing other hypotheses not 
related to curvature. This seems to be a challenging problem; cf. [1] for some further discussion. 
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